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The two-dimensional (2D) continuous elastic energy model for isotropic tubes is reduced to a 
one-dimensional (ID) curvature elastic energy model strictly. This ID model is in accordance with 
the Kirchhoff elastic rod theory. Neglecting the in-plane strain energy in this model, it is suitable 
', to investigate the nature features of carbon nanotubes (CNTs) with large deformations and can 

I ' also reduce to the string model in [Z.C. Ou-Yang et al, Phys. Rev. Lett. 76 4055 (1997)] when 

' the deformation is small enough. For straight chiral shapes, this general model indicates that the 

, difference of the chiral angle between two equilibrium states is about it/6, which is consistent with the 

lattice model. It also reveals that the helical shape has lower energy for per atom than the straight 
shape has in the same condition. By solving the corresponding equilibrium shape equations, the 
helical tube solution is in good agreement with the experimental result, and super helical shapes 
are obtained and we hope they can be found in future experiments. 



, . I. INTRODUCTION 

I ' In the past two decades, the CNTs initially synthesized by lijima attracted many researchers' attention due to 
^ . the excellent physical characteristics and potential applications in many apparatus and nano-instruments, such as field 
emission sources 0], probe tips 0,3 and quantum wires Q. Experiments indicate that the configurations of CNTs, 
such as their radii, lengthes and helicity, strictly determine their physical capability, but it is difficult to precisely 
' control those configurations during production processes. So, although the CNT is a particularly important functional 
nanomaterial, obtaining macroscopical bulk materials is a challenging problem Q. Recently, an important progress 
reported by Davis et al. shows a way to obtain macroscopical fibres of the SWNTs using the self-assemble method 
in chlorosulphonic acid. Moreover, the mechanical parameters of CNTs are not unification due to the size effect, 
which confines their applications as a reliable high strength material. For instance, the Young's modulus will decline 
O with the increase of the diameters of CNTs @ . The above problems indicate that there are still some challenges to 
1^1 , overcome before we can widely put CNTs in practice. 

An interesting phenomenon is that CNTs often present as beelines (including zigzag and armchair shapes (9l-[ll|). 
in helixes [H, M, fisj as well as rings [l^, ■ And a possible reason for those shapes is the thermodynamic effects in 
^ ' different synthesizing methods. Particularly, periodic defects (heptagon and pentagon cells) play an important rule in 
the forming of helical shapes [H, [HI • SWNTs are generally taken as a frizzy graphite layer and multi- walled carbon 
nanotubes (MWNTs) consist of multiple rolled layers of graphite. Many physical properties of the CNT are obtained 
by calculating the interaction between its carbon lattices. There are also many researchers who take the CNT as 2D 
• , continuous tubules [T9l - [2^ . Ou-Yang et al. jl^ continued the lattice model which provided by Lenosky et al. [1^ 
' and pointed out that SWNT's free energy is similar to the model for vesicles [l^] when ignored the in-plane strain 
energy. In Refs. [21, . ,22.] Tu and Ou-Yang provided a general 2D model which considered the in-plane strain energy 
and revealed that the effective Young's modulus of MWNTs dependents on the layer number. 

Although the CNT is taken as the ID elastic material, by far there is still not a strict model to connect the lattices 
model and the Kirchhoff elastic rod theory. The string model |1S] provided a recommendable way to make up this 
missing link, but it cannot be used to investigate the mechanical behaviors when outside forces act on CNTs, because 
^ , the in-plane strain energy can not be ignored in this circumstance. So it needs to construct a complete ID model 
^ _' which should contain the in-plane strain energy to connect the lattices model and the Kirchhoff theory. Besides, 
although the lattice model can geometrically tell us that the defects induce the helical shapes, it lacks a reasoned 
physical theory to expound the reason of the emergence of helical CNTs. Can the continuous model explain why 
there are so many helical CNTs? Further, considering that ID structures can form super helical shapes, such as DNA 
chains, so can wc obtain super helical CNTs? These questions need to be investigated deeply. Not only will they 
help us to understand the physical characters of the low dimensional systems, they will provide us with new materials 
and methods to design nano-instruments. In this paper, we will give a complete ID continuous CNTs model which 
contains the in-plane strain energy to discuss the above problems. This paper is organized as follows: In Sec. |lTl the 
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2D continuous elastic shell model is reduced to a ID curvature elastic model strictly. This complete ID model is in 
accordance with the isotropic Kirchhoff elastic rod model and is suitable to investigate the mechanical behaviors of 
CNTs. In Sec. Illli a concise model which ignores the in-plane strain energy is used to study the nature features of 
CNTs. By solving the corresponding shape equations, the helical solution is in good agreement with the experimental 
result and super helical shapes are obtained. Finally, a short discussion is presented in Sec. HV] 



II. A COMPLETE ID MODEL FOR ISOTROPIC ELASTIC TUBES AND CNTS 

Let the central line of an elastica tube be R = R(s), c? = R be the tangent vector (an overdot denotes differential 
with respect to s which is the arclength of the central line of the tube), f3 — 11/ K be the main normal vector 
and 7 = Q? X /3 be the binormal vector, between those unit vectors, there are the Frenet formulaes: da/ds = K/3, 
d(3/ds — —Ka + and dj/ds — ~Tj3, where K and r are the curvature and torsion [25j of the central line R, 
respectively. The shape of the tube can be obtained by this way: the central point of a ring with the radius rg moves 
along a line R and keeps the ring upright to the tangent of R. It means the ring is in the normal plane of the central 
line R. Let Y be the shape of the tube, there is 

Y = R + r = R + ro cos 6'/3 + ro sin 6*7, (1) 

where the parameter Q < 9 < 7r/2. Using the Frenet formulaes, we obtain the mean curvature and Gaussian curvature: 
H ~ „ ^7"^"^?/^ '^°a^ t\ and A = — fc cos e respectively. Here we give a method in appendix!^ to show how to attain 

2ro (ro-K cos w — 1) ro(ro-K cos p — 1) ' ^ ^ ' — ' 

these two curvatures. The 2D curvature elastic energy of a tube is [2l|, [26l - [29l | 



24(1 - 1/2) 



^ [{2H)^ - 2(1 - i^)A] da, (2) 



where Y is the Youngs modulus, v is the Poisson ratio, h is the thickness of the tube and da is the area element. The 
above energy derives from bending the graphene to SWNTs. For any compact, closed 2D surface, / Ada is a constant 
because there is the Gauss-Bonnet theorem (Sli 



Ada = 2TrEr, (3) 

where Er is the Euler characteristic which only depends on the topological structure of the surface. For spherical 
topology surface, there is Ej. = 2. For torus, cylinders and tubes with infinite length, there is Er = 0. In the following 
text, we will study the torus and tubes with infinite length, so we eliminate the constant term associated with f Ada 
in ©. Considering da = ro(l — roK cosd)dd and r^K^ < 1, and using the Euler integral, Eq. 1^ is reduced to 

= n ■.2^ / o d'- (4) 



12ro(l- 1/2)7 ^1 - r^K^ 
The in-plane strain energy can be expressed as (26l-[2 



Yh 

^» = :J\ ^ [(^- + ^y)' " 2(1 - iy)(e.,.ey - ely)] da. 



2(1 



xyj 

(5) 



where Sx, Sy and Sxy are the axial, circumferential, and shear strains, respectively. In Fig. (TJa), we show a fragment 
of SWNT with the radius ro and the length of its central line is ds. We suppose that its central line has a bend 
angle d(t), so the radius of curvature is R = \/K — ds/d(f>. Taking the SWNT as an isotropic elastic tube and the 
central line is nonretractable, and supposing that the initial state is a straight tube, from Fig. [T]we can see that the 
protraction of the point P along the axis direction is dso — [R — ro cos{9 + d9)]d(l) — ds ~ —ro cos 0d<j). Then the tensile 
strain on point P is e^: = ^ — —roK cosO. We simply choose the circumferential strain Sy = —ve^ = vr^KcosQ. 
Correspondingly, the shear strain on point P is £xy — ^''o^- Note ^ = + there is Exy = tq(t -|- x)/2. Here the 
definitions of d and x please see the caption of Fig. [T] Using the above results, Eq. ([SJ is reduced to 

■nrlYh f p 2 , / , -^21 



F. - 2^^^ / [(1 + iy)K' + (r + x)'] ds. (6) 
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FIG. 1: (a) A fragment of a curving SWNT with the radius ro, the length of central line ds and the bending angle dtf), C and 
C" are the cross section rings on the two ends, (b) The left cross section of the SWNT. (3 and 7 are the main normal vector 
and the binormal vector of the central line, respectively, [x, y, z) are the local coordinates with the x, y axes fixed on the cross 
section ring and with the 2 axis superposed to the tangent of the central line, x is the angle between (3 and the x axis. Point 
P is fixed on the cross section ring (on the left end) and the angle between OP and /? is 9. When the cross section ring moves 
along the central line from the left end to the right end, P changes to P' (Note that the cross section ring can turn around the 
central axis.). In this process, the angular displacement of point P in the local coordinates (x, y, z) is dO and the change of x is 
dx- Supposing that the initial state is a straight tube, this two charts indicate the protraction of the point P along the central 
line is dso = —rocos9d4>, the corresponding tensile strain is Ex ~ ^ — —roK cos6 and the shear strain is Sxy — ro{r + x)/2. 



Then, the ID energy density of a tube can be written as 



nYh^ 1 Trr^yh 2 

12ro(l - 1^2) yT=7p^ + 



2(1 + 

For a multilayered tube, the energy density is 



(r + x)- (7) 



r {F/b + 2nrag/b)dro, (8) 
•I Pi 

where b is the distance of two neighbor layers, g is the surface energy density between two neighbor layers, pi and po 
are the inmost and outmost radii, respectively. Specially, when b = h, there is 



Po{l + ^l-K^pi 



]^YlrK^ + ]^GIt{T + xf + PS^, (9) 

where 1^ — T^ipt — Pi)/4: and It = T^ipg — Pi)/2 are the moments of inertia of the cross section around its diameter 
and central axis, respectively. G — Y/{2 + 2v) is the shear modulus, P = g/b can be taken as the volumetric energy 
density and Sq — Tr{p1 — pf) is the area of the cross section of the tube. When pi = 0, the second and the third terms 
in the right hand side of the above equation compose the typical isotropic KirchhofF elastic rod model. 

In Refs. [2I, 22], Tu and Ou-Yang continued the lattice model and showed that the energy density for SWNTs is 
similar to the elastic shell model. Using their results and considering Eq. (O, the energy density for SWNTs is 

where A = 7rfcc/ro, B = {I ~ iy'^)nr^kd, C = B/{1 + v), kc = 1.62 eV and kd = 22.97 eV/A^ [22^. Considering 
rQK^ <C 1, r-Q = 0.5 ~ 10 nm and B ^ TrkcVo, the first term on the right hand side of the above equations can be 
ignored and similar models are discussed in Refs. [s^, HI] to study the mechanical properties of nanosprings. The 
energy density for MWNTs is 



J^mw ^ / iJ^sw/b + 2T:rog/b)dro, (11) 



FIG. 2: A curving SWNT is composed of several straight segments connected by abrupt corners. The average length for each 
segment is AS and the average corner between two neighbor segments' central lines is A<j!>, so the curvature for the central line 
is K~ A^/AS. 



where b ~ 0.34 nm and g = —2.04 cV/nm^ [33]. Specially, for a straight double- walled carbon nanotube (DWNT), 
there is = t = in (fTO|). Moreover, we should note that x is just the chiral angle in this case. Considering each 
straight layer often has fixed chiral angle, there is x = (actually, this result means there is no strain energy, similar 
result will be shown in the later text), thus the energy density for a DWNT can be written as 

-^d^ =^fc, f — + ) +27r(ro + 6/2).9, (12) 

where Tq is the inner radius. Choosing kc = 1.4 eV [s^, J-dw = yields 

ro = 6.8 A. (13) 

For dense SWNT ropes, supposing that each SWNT is continuously enveloped by its neighbors, the model in (fT2)) 
can give the optimal radius for SWNTs and the similar result obtained by Zhang et al. is 6.8 A |35l |. 

In the above discussions, the energy density in (jlOp which contains the contribution of the in-plane deformation is 
suitable to study the mechanical behaviors of SWNTs. However, if we want to study the natural features of SWNTs 
without any outside forces, the in-plane deformation terms in (|10[) should be neglected. The corresponding model will 
be studied in the following section. 



III. ID MODEL FOR CNTS WITHOUT IN-PLANE STRAIN ENERGY 



A. Analytical results 



According to the experimental observation [l2|, |36| , helical SWNTs with periodic defects can be treated as the shape 
in Fig. [21 Seeing Fig. [21 a SWNT is composed of several straight segments connected by abrupt corners. The average 
length for each segment is AS* and the average corner between two neighbor segments' central lines is A^, so the 
curvature for the central line is -fC ~ Ac/)/ AS. Supposing that the average torsion angle around the a axis between 
two neighbor segments is Ad ( A9 here is similar to d9 in Fig.llJ, there is r ~ {Ad — Ax)/ AS. However, although the 
central line of the shape in Fig. [21 has curvature and torsion, which are due to the abrupt corners not to the in-plane 
strain. Thus the second and third terms on the right side hand of Eq. (|TOl) should be ignored for SWNTs with out 
outside forces. Seeing an example, for a helical shape with K = j^i^f^i and r ~ B!^%h'^ ; define lo = h/R, there is 

A(t)^u}{A9- Ax), (14) 

here Ax is the chiral angel difference between two neighbor straight segments. In [l2|, Zhang et al. found Atj) — tt/6, 
and particularly in Fig. 5a of [13], we can see that two neighbor segments are composed by a zigzag shape and an 
armchair shape, which indicates Ax = — 7r/6. Moreover, Ou-Yang et al.^^ pointed the shape in Fig. 2a of [12] 
satisfies a; = 1, then substituting the above results into Eq. (ITl)) yields 

Ae = 0. (15) 

This result indicates t + x— A0/AS — in Eq. (jTU]) strictly. The length for per helix turn is Sp — I-k^J + h? , and 
the segments number for per helix turn is 
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For = 7r/6 and u! = 1, there is ~ 9. This result is close to the experimental observation that there are about a 
dozen bends per helix turn ^12] . 

According to the above analysis, to investigate the holistic nature features of SWNTs, the in-plane strain energy 
should be neglected. Thus the total energy density for SWNTs can be written as 

where A is the Lagrange multiplier. As to SWNTs, A can be taken as the average intensity of the effect between 
a SWNT and it's neighbor, such as the dense SWNT ropes model in [s^. For MWNTs, we can take A as the line 
tension coefficient due to the effect between one layer of MWNTs and its two neighbor layers. Thus the above energy 
density in ([TTl) also is suitable to each layers of MWNTs. Comparing Eqs. ([TT|) with (fTT]) . we have A ^ ^irrog ^ —10 
eV/nm and ^jp- ^ — 1. This result will be proved by our later calculations. 

One can easily find that the model in ([TT]) will reduce to the string model when r^K^ <C 1. Moreover, from ([2]) 
to dl]), we need r'^K^ < 1 not r^K^ <^ 1. This small difference will make the model in ([T7| is suitable for the CNTs 
shapes with lager deformations and give us much more abundant shapes than the string model. For a ID elastica 
structure with the energy density functional J- = J-{K,t), the equilibrium shape equations have been obtained in 
[37I [38| by discussing the first variation of the energy 5{§ Fds) = 0. Making use of their results, we attain the 
equilibrium shape equations for the energy density in p7|) 



rl (1 + rlK"^ - 2r^K^)k + Sr^/C (3 + 2rlK'^)k'^ 
^K{l-rlK^f{l-2rlK^+rlT^) 

-\K{l~rlKy^ ^Q, (18) 
2kT{l + 2rlK'^)+Ki{l-rlK'^) = Q, (19) 

where we define A = A ring solution with the radius R yields 

(1 - 2r2/i?2) _ A(l - r2/i?2)3/2 ^ q. (20) 

According to the experimental results, in most cases there are 0.5 nm < < 5 nm and 100 nm < R < 500 nm 
[H, [l3|, so we have r^/R^ — >■ 0, which means A ^ 1 in Eq. (l20l) . Specially, when A = in Eq. ([20)) . we obtain 

R = V2ro. (21) 

This Clifford torus solution for vesicle was found by Ou-Yang (39| and proved by a coming experiment (ioj . Avron 



and Berger [4l| gave some details about the torus nearby -R/rp = 2. So the Clifford torus with R/r^ — \f2 is easy to 



be constructed, such as the shapes in series II of Fig. 3 of [4^, which are close to this shape. 

Substituting the helix solution K = Kq = -^2^^,t = tq = ^^2^112 into Eq. and defining oj — h/R, i] = r^/R 
and the helical angle tj^h = argtanw, we have 

+A[(l + c.2)2-^2]3/2^Q^ ^22) 

We show an example solution in Fig. [31 which is consistent with the values in Fig. 3(b) of [l^. Specially, the zero 
energy state w = 1 (V'h = 45°) [H] yields A = -2/^/4~^. If A = 0, Eq. ([22]) is reduced to 

(23) 



Considering < 77 < 1, this equation indicates ^ < 77 < 1, < {^/T3 — 3)/2 ~ 0.3 and the helix angle tph < 28.8°. 

However, when ^ < i] < 0.764, the shapes are self-intersected. The valid region is 0.764 < 77 < 1 (14.1° < < 28.8°) 
and these shapes are close to the Cioso shape in Fig. 1(b) of (il j. 

Particularly, if 77 = 1 and w — ;> 00 in p2p. we obtain a cylinder solution. However, if 77 = 1 but uj is finite, what 
kind of shape we can obtain? If so, we will see that it is nothing but the chiral configuration. Supposing that the 
state with 77 = 1 and — 00 is the zigzag shape which has two sp^ bonds of each carbon hexagon paralleled to the 
axial line of the SWNT and has the chiral angle V'/i = 7r/2, so the armchair state should be with = 7r/3 (Note 
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that the chiral angle tpi^ in this paper is the complement angle for the definiens in [45|). Here we define the reduced 

ro_ ttO 



total energy density for equilibrium helical shapes: ft — -^F^ , using the method in j43|, there is 



(l-w4)r;2 



[(l+w2)2-7y2]3/2' 



(24) 



where J^i — = and J^" = J^|/i:=ifo_T=To- For straight shapes with r/ = 1, equilibrium condition dfl/dio — 

yields 



io^oo, Lu^ Y ^(10^^^ + 10^^^ + (25) 

These correspond to 4'h = 7r/2 and i/jh = 58.2° sa 7r/3, respectively. We show il = il.{tph) in Fig. |4]which indicates 
that the chiral angle difference between two equilibrium states is about 7r/6. In fact, the 2D continuous model in 
pT)) can not tell us which state in ((25)) is the zigzag shape. However, it gives us the angle difference between two 
equilibrium states. This angle difference Atpi^ = 31.8° ~ tt/6 is in good agreement with the lattice model. Further, 
considering that each carbon hexagon has the area sq — 3\/3tP /2 and possesses two carbon atoms, where d — 1.42 A 
is the equilibrium distance of two neighbor carbon atoms, the energy for each carbon atom is 

E^^s^^'^^'f , (26) 

47rro 8rl ^(1 + cj2)2 _ ^2 ' 

where the unit for tq is A and A = 0. Then choosing = 1.62 eV, [ll we have E = and E = 2- 20869 

equilibrium states in (note = 1), respectively. We simulated the corresponding values in 1461] obtained by QMD 
method and found that the energy for each carbon atoms satisfies E = 2-°^J^^ for zigzag shapes and E = for 

armchair shapes. Clearly, the result obtained by the continuous method in (|26p is close to the corresponding result 
obtained by QMD method, and it also indicates the two states in (|25p are the zigzag shape and the armchair shape, 
respectively. If we want that the result in ([^51) can consistent with the QMD result in 46] more exactly, we can choose 
kc = 1.55 eV. Moreover, for a chiral shape with rj = 1 and cj 7^ 0, and a helical shape with 77 < 1 and the same oj 
as the chiral shape, Eq. (^5)) indicates the helical shape will have lower energy for each atom. It explains why there 
are so many helical shapes are found in experiments. The above results indicate that the continuous model in (jl7p 
which based on taking the sp^ bonds as the geodesic lines on a tube also contains a little information about the lattice 
structure. 



Numerical results 



Let the central line of the SWNT be R = {Rx, Ry, Rz}, Rx = cos cos (^9, Ry = cos sin 93 and Rz = sin6', where 
= 6{s) and (p = ^p{s) are two Euler angles with the variable s, there are 



i^2 ^ ^2 ^ ^2 ^,Qg2 

TK^ = ^(2^2 + ^2 ^Qg2 Q-^ gjjj Q (^^Q _ Q^^ Q 



(27) 
(28) 
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FIG. 4: The chart of the reduced total energy Q, = Q{iph) in (|24l) with r] — 1 for chiral shapes. It indicates that the chiral angle 
difference between two equilibrium states is Atph = 31.8° ~ 7r/6, which is consistent with the lattice model. 




FIG. 5: A positive super helical shape with the initial conditions: ro = 1 nm, A = —0.978, 6'(0) = ip{0) = ifi{0) = 0, ^(0) = 0.04 
nm "\ ip{0) = 0.26 nm"^ and ^(0) = 0.01 nm . The first helix is right-handed and the second helix is left-handed. 



Substituting the above expressions into Eqs. (|T8|) and (|T9|) . we obtain two tedious third order equations about 9 and 
if (see appendix IbJ . Solving this two equations, we obtain several interesting shapes. Fig. [5] shows a positive super 
helical tube, Fig. [5] shows a negative super helix tube and Fig. [7] shows a right handed helical ring. Experimental 
basketwork in Fig. 2(b) of [47[ which contains many super helical MWNT shapes has excellent mechanical capabilities. 
As to these unattached super helical shapes in Fig. [5] and Fig.[6l we can conclude that they have strong and restorable 
retractility like the DNA chain. So, they are perfect functional materials and have large potential in constructing 
nano- instruments. 



FIG. 6: A negative super helical shape with the initial conditions: ro = 3 nm, A = —0.97, 6(0) = (^(0) — 0, 9{0) = 0.018 nm 
ip{0) = 0.08 nm"\ 9(0) = 0.0009 nm~^ and if{0) = 0.0011 nm~^. The first and the second helixes are right-handed. 
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FIG. 7: A right-handed hehcal ring with the initial conditions: ro = 3 nm, A = -0.96, 61(0) = ^3(0) = 0, 9(0) = 0.002 nm 
ifiiO) = 0.07 nm-\ 8(0) = 0.002 nm"^ and ifi{0) = 0.00261 nni'^ 

In the planar case, the valid shape equation is reduced to 



where T = tqK. This is a non- linear equation and it is difficult to be solved generally. We numerically solved this 
equation but all the solutions we obtained are similar to the shapes in Fig. 4 of [1^. 



In conclusion, we have shown a connection between the elastic shell model and the Kirchhoff elastic rod model. 
Combining Refs. [l^[lH and this work, a complete method to deal with CNTs has been constructed. From lattices 
model to 2D continues elastic shell model and further to ID continues elastic rod model, this method gives a recom- 
mendable approach to del with the multi-scale low dimensional systems. The ID model in ^ which contains the 
in-plane strain energy is suitable to investigate the mechanical behaviors of CNTs. But we should note that the first 
term on the right hand side of ([9]) is the distinct difference, when compared with the isotropic Kirchhoff elastic rod 
model. This difference will make CNTs have unusual mechanical behaviors which have not been known. Moreover, 
we should note that the model in (jTU]) is obtained by taken CNTs as isotropic perfect tubes. If there are plenty of 
pentagonal and heptagonal defects for bended CNTs, this model seams need to be changed. A simple way to adapt 
this change is to adjust the constants kc and kd, such as the work in [33] where the authors chose kc — 1.4 eV which 
is different to the value kc = 1.62 eV in p^ . 

For straight chiral CNTs, our study indicates that the difference of the chiral angle between two equilibrium states 
is about 7r/6, which is consistent with the angle difference between the zigzag shape and the armchair shape obtained 
by lattice model. Our study also reveals that, if a helical shape and a straight chiral shape have the same radius rg 
and the same cj, the former structure will have lower energy for per atom than the later one, which explains why there 
are so many helical CNTs in experiments. Since there are super helical solutions for the equilibrium shape equations, 
we hope they can be found in future experiments and the super retractility also can be found for these shapes (If 
it is only the string model fl^, there are only helical solutions but not super helical shapes [49]). How to produce 
super helical CNTs in experiments? Yin et al. [50| provided a way to construct the super CNTs using a cylindrical 
template. If we choose a helical template, super helical CNTs may be available. 

Finally, we would like to point out that how to use the lattice model to construct the super helical CNTs needs 
further discussion, which will be our future work. 



rg(l + - 2T^)f + 3r^T(3 + 2T^)t 




(29) 



IV. CONCLUSIONS 
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Appendix A 

Defining (/i, v) to be the local coordinates on the tube Y, where is along the tangent of the central line of the 
tube R and v is upright to R, there are dfi = ds and di> — r^dQ. Making use of the Frenet formulae, we obtain 

= Y^-Y^ = (1- roi^ cos 0)2+ r^T^, (Al) 
F = Y^ • Y, - ror, (A2) 
G = Y, • Y„ = 1, (A3) 



where Y^; = dY jd\x. The main normal vector of the tube is 



X Y„ 
|Y^ X Y, 



= r/ro. (A4) 



Consequently we obtain 



L = Yf,f,-n::^ Kcose{l-roKcosey -roT^, (A5) 
M = Y^„ 11 = -T, (A6) 
N = Y,,-n = -l/ro. (A7) 



The mean curvature and Gaussian curvature are 



LG-2MF + NE _ l-2roKcose 
^ ~ 2{EG-F^) " 2roiroKcos9-l)' ^ 
^ _ LN-AP _ Kcose 

EG-F^ ^ ro{roKcos9-l)' ^ ' 



Appendix B 



In this part, we show the shape equations (|18p and (|19p with the Euler angles as variables. Substituting expressions 
(P7)l and (|^ into Eqs. ([T^ and we obtain two tedious equations 



4:rl(6ip'^ sin 6* cos 6* - 
+4rl { 1 + (^2 + ^Qg2 61) [l - 2ro 



- ipip cos t/j 

X + (^2 ^Qg2 g,^|-gi2 _^ ^2 ^^^g2 0_q2^2 ^9 

-9tp'^ sin 61 COS 61 - 29tp^ sin 26* + 99^^^ + tptp^^^ cos^ 

-4 [1 - r2 (02 + ^2 C0S2 0)] ' X { (02 + ,^2 ^^^2 g)2 
1 - 2r2(02 + ^2 ^Qg2 0) + A(l - r2(02 + ,^2 ^^g2 g)^3/2- 

+ro [202 sin + V? cos 0(^32 sin cos + 0) - 0(/3 cos 0] ^ | 

= 0, 



(Bl) 
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{96 + ipifi cos^ e - Oifi^ sin cos 0) [l + 2rl [d^ + (p^ cos^ 
X 2 [2^^(/3 sin 6* + V3 cos (?((^^ sin 6* cos 6* + Q) - 9(p cos 9\ 

+ [1 - + (p^ cos^ 61)] 1 cos 6» + 3^^<^ sin 6i 
+0^ cos^ 6{'f'^(p sin 6 cos 6- 2e<p + <pe^^^) 
cose{2e(p + ipe^^'^) + [0^(11 cos6i + cos 36i) 
-Ae^p sin 61 - 4^.(3) cos 6»] + </9^ cos e [3</9^6» sin 60056- 2P 
+ cos^ 6'(v)* cos^ 6 + 2ip^ - I 

= 0. (B2) 

Solving the above equations, we get = 0{s) and ip = </?(s), and consequently obtain = cos 9 cos ipds, Ry = 

s s 
Jq cos6svcnpds and Rz = Jg sin^ds. The SWNT shape can be written as 

Y = {x,y,z}, 

X = {(P' + ijp^ cos^ 6) ~ [ro sin Lp(d sin (j) — ip cos cos 

— ro sin^^ cos cos (/) + cos sin (/))] +-R2:, (B3) 
F = (^^ + cos^ 6*) ~ [ro cos Lp{^p cos 9 cos (/) — ^ sin 0) 

— rosin^sin(p(^cos(/) + (yicos^sin^)] + Ry, (B4) 
Z = ro cos 9{6'^ + </j^ cos^ 6*) ~^ {9 cos (f) + ip cos 9 sin (/)) 

+ Rz- (B5) 



